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ABSTRACT 

As planets form they tidally interact with their natal disks. Though the tidal perturbation induced 
by Earth and super-Earth mass planets is generally too weak to significantly modify the structure 
of the disk, the interaction is potentially strong enough to cause the planets to undergo rapid type I 
migration. This physical process may provide a source of short-period super-Earths, though it may 
also pose a challenge to the emergence and retention of cores on long-period orbits with sufficient 
mass to evolve into gas giants. Previous numerical simulations have shown that the type I migration 
rate sensitively depends upon the circumstellar disk's properties, particularly the temperature and 
surface density gradients. Here, we derive these structure parameters for 1) a self-consistent viscous- 
disk model based on a constant a prescription, 2) an irradiated disk model that takes into account 
heating due to the absorption of stellar photons, and 3) a layered-accretion disk model with variable 
a parameter. We show that in the inner viscously-heated regions of typical protostellar disks, the 
horseshoe and corotation torques of super-Earths can exceed their differential Lindblad torque and 
cause them to undergo outward migration. However, the temperature profile due to passive stellar 
irradiation causes type I migration to be inwards throughout much of the disk. For disks in which 
there is outwards migration, we show that location and the mass range of the "planet traps" depends 
on some uncertain assumptions adopted for these disk models. Competing physical effects may lead 
to dispersion in super-Earths' mass-period distribution. 
Subject headings: Planetary Systems: Protoplanetary Disks 



1. INTRODUCTION 

New high-precision radial velocity and transit sur- 
veys allow observers to detect low-mass and modest- 
size planets, particularly when those planets with short- 
period orbits. Evidence is emerging that these close-in 
super-Earth mass planets are significantly more ubiqui- 
tous than the ir hig her mass cousi ns, the hot Jupiters 
([Mayor et al.l 120091: iMavor et all 120111 : iHoward et all 
2010). 

Most theoretical models of hot Jupiter formation as- 
sume that planets form far from their parent stars 
and migrate inwards , either by planet-planet scattering 
dRasio fe Fordl 119961: iWeidenschilling fe MarzaTH 119961: 
ILin k Idalll997HMurrav et al.lll998ft. secular interactions 
such as the Kozai resonance (IWu fe Murravll2003D, sec- 
ular chaos (iNaeasawa et al.1 120081: iMarchi et al.l [2001 
IWu fc Lithwickll2011[ ), or thr ough tidal intera ctions with 
a gaseous circumstellar disk (|Lin et al.l I1996T) . This mi- 
gration is required for hot Jupiters because it is difficult 
for them to form in situ. The gas disk, which must be 
present to supply material to these planets, damps the 
eccentricity of planctesimals and embryos. This means 
that the growth of a protoplanetary embryo, on a cir- 
cular orbit, stalls once it acquires an isolation mass 
(Mj SO ), i.e. it has captured all solids within a feeding 
zone centered around its orbital radius r. The width 
of the feeding zone is approximately ~ lOrnm where 
the Hill's radius for a planet of mass M p around a star 

of mass M* is r H iii = r[M p /(3M*)] 1 / 3 (|Kokubo fc Idal 
l2002f) . This leads to an isolation mass in a disk with the 
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local surface density of dust and planetesimals (E^) is 
M iso ~ (E d a 2 /M*) 3 / 2 M*. 

If we adopt a value for E^ inferred from the empir ical 
minimum mass nebula model (MMN, lHav ashi 198U see 
below), Mj so would be a fraction of Martian m ass at 1AU 
and a few times Earth mass (M e ) at 5AU (llda fc Lid 
2004). This maximum mass increases with the distance 
from the ce ntral star as th ere are more solids available 
at large r (jLissauerl Il987t ) . In the conventional core- 
accretion scenario, efficient gas accretion onto gas gi- 
ants needs to be preceded by the emergence of core s 
with masses in excess of - lOMq JPollack et al J 11996ft . 
Progenitor cores with such large masses can only be as- 
sembled beyond the snow line in an MMN. 

The same argument also suggests that it is difficult 
to form the rich population of known close-in (within 
~ 0.1AU) super-Earths (a few to 15 M®) in situ. In order 
for there to be an adequate amount of building block ma- 
terial to assemble these planets the magnitude of E^ must 
be more than an order of magnitude larger than that 
of the MMN. If this necessary condition is marginally 
satisfied then, in the presence of gas, instead of produc- 
ing a single larger embryo, many embryos with much 
smaller isolation mass would emerge. Under these con- 
ditions some rocky Earth analogs or super-Earths could 
have formed in situ through collisional mergers after the 
gas disperses (as it can no longer damp embryos' mutu- 
ally excited eccentricities). But giant impacts could still 
cause a substantial loss of volatiles in post-gas evolution, 
so this method may not lead to the emergence of the 
known low-density super-Earths. 

In order for the magnitude of M- mo to reach a few M ffi 
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at 0.1 AU, the local £d must be increased by two orders 
of magnitude over the MMN. In disks with a solar com- 
position, the corresponding gas surface density (£) would 
be sufficiently large that gravitational instability would 
induce a rapid redistribution of mass and angular mo- 
ment um (jDurisen et al.ll2007t[Gammiel[200ltlBolev et"aH 
120061) . Protoplanetary embryo formation in such disk re- 
gions may be challenged because the disk temperature 
would be above the condensation temperature of most 
refractory grains and the £ would decline substantially 
before the planetesimals can grow significantly. In gravi- 
tationally stable disks with a £ distribution comparable 
to that of the MMN, the dust to gas ratio would need to 
be enhanced by mo re than an order of magnitude relative 
to its solar value (|Ida fe Linll2005D . This enhancement 
can be induced by planetesim al drift (e.g. iKretke et al.1 
120091: lHansen fe Murravll20TlT ) or by embryo migration. 

Planetary embryos migrate due to their interaction 
with gas in their natal disks. Their tidal perturbation 
excites density waves in their Lindblad and co-rotation 
(or horseshoe) resonance regions. These waves carry an- 
gular momentum flux which is determined by the sur- 
face density (£) and tem perature (T) of the disk ga s 
at the excitation location (jGoldreich fc TremaineilT980D . 
These waves dissipate and deposit angular momentum 
to the disk as they prop agate away from the resonances 
(|Papaloizou fc Linl I1984D . The disk responds with an 
evolving £ a nd T distribution wh ich introduces a feed- 
back effect (jTakeuchi et al.l I1996D . In order to calcu- 
late the full non-linear disk response, it is necessary to 
perform numerical hydrodynamic simulations of planets 
embedded in two or three dimensional disks. And mu- 
tual interactions between multiple planets embedded in 
a disk can lead to a wide varieties of diverse behaviors in- 
cluding outwards migration and migration stalling (e.g . 
IMasset fc Snellgrovdl200l1; IPodlewska-Gaca et al.ll2012ti . 
Nevertheless, if we restrict our investigation to planets 
with modest masses, the disk structure is not strongly 
perturbed and their net linear Lindblad and co-rotation 
torque can be calculated by summing up the contribution 
of each distinctive resonance. In disks with a power-law 
£ and T profiles it is possible to describe the net torque 
on the planets as a linear combination of the radial £ 
and T gradient (jTanaka et alll200"2t iPaardekooper et ah! 
I2010al) . 

In this paper, we are primarily interested in super 
Earths' type I migration. The mass of these planets 
is generally too small to significantly perturb the struc- 
ture of the disk. In such limit, super Earths' migra- 
tion rate can be deduced from the structure of unper- 
turbed disk. Here, we apply various disk models to 
some recently obtained prescrip tions for type I migra- 
tion (jPaardekooper et all l2010bl hereafer PBK10) and 
determine the direction and rate of migration for planets 
with a range of masses. 

The simplest empirical prescription for the disk struc- 
ture is the MMN model which is based the assump- 
tion that all planets in the solar system attained their 
mass in situ and they accreted all the heavy elements 
in th eir neighborhood ( Wcidcn schillinel Il977bt IHavashil 
119811) . Through the augmentation of volatile elements for 
a solar elemental distribution, we can infer a £ oc r~ 3 / 2 
distribution. With an additional assumption that pro- 



tostcllar disks are heated to their local equilibrium tem- 
perature, the disk temperature T oc r -1 / 2 . The MMN is 
a reasonable fiducial model, albeit the in situ formation 
assumption is inconsistent with the current concept of 
planetary mobility. 

Alternatively, we can construct a set of dynamical 
models for the protostcllar disks in which the £ and 
T distribution are determined with physical principles. 
Due to uncertainties in the mechanism of angular mo- 
mentum in disk it is commo n to use the a- viscosity model 
(jShakura fc Sunvaevl I1973D in which one assumes that 
the disk evolves due to a viscosity proportional to the 
sound speed (c s ) and disk scale height (/i), 

v = ac s h. (1) 

For the disk models we adopt here we first present re- 
sults for disks with the common assumption that a is 
constant. We then investigate disks in which we assume 
the magneto-rotational instability (MRI) excites turbu- 
lence and leads to angular momentum tran sfer in regions 
with sufficiently high ionization fraction (jHawlev et all 
I1995D . In most planet forming parts of the disk (ranging 
from a fraction to a few 10 or AUs), MRI active regions 
are confined to the surface layers which are exposed t o 
ionizing stellar photons and cosmic rays (|Gammi 
In order to include this effect, we generalize the quasi 
steady state accretion disk model to allow for a variable 
a parameter. 

There are a number of groups who have presented 
sophisticated str u cture models of MRI- active disk s 
(ISano et all l2000t iTeraueml [20081: IKretke fc LiiJ l2010ft . 
Due to the uncertainties on both MHD turbulence and 
dust grains' size distribution in the disk surface layer, we 
choose to use a model with fewer assumptions to simplify 
this issue. 

A few studies have already begun investigating planet 
migration in MRI-active disks includ ing t he effects of a 
dead zo ne. iKuchner fc Lecail ([2002D and IMasset et all 
(|2006bl) investig ated planet stallin g at the inne r edge of 
the dead zone. IMatsumura et all (|2007l I2009D studied 
the growth and migration of Jovian planets in disks with 
dead zones. All of these studies indicate that planet mi- 
gration is sensitive to the structure of the disk. In this 
paper we present the evolution in simple disk models in 
order to discuss the generic disk profiles that can cre- 
ate "planet traps" and how these are likely to manifest 
themselves in physical disks. 

As we have indicated above, the rate and direction of 
type I migration is determined by the competing con- 
tribution from differential Lindblad and corotatation (or 
horseshoe) torque. In our computation of the condition 
for corotation torque saturation, we are mindful that the 
horseshoe region has a toroidal shape. Its height normal 
to the disk plane is comparable to its width in radial di- 
rection. For planets with mass less than a few M©, the 
entire horseshoe region may be embedded in the dead 
zone. It is not clear whether there is any significant vis- 
cous flow across the horseshoe zone to prevent satura- 
tion. In these low-viscosity regions, embedded planets' 
long-term interaction may also modify the £ distribu- 
tion, tri gger secondary instabilities, and halt their mi- 
gration iBalmforth fc Korvcanskvl 120011: iLi et all I2005L 
120091: lYu et alll2010D . 
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This paper is organized as follows. In wc briefly 
recapitulate relevant migration prescriptions. In ££3j we 
construct three sets of disk models which take into ac- 
count surface irradiation and the presence of the "dead" 
midplane between MRI active surface layers. In §2J we 
present the results on how planets will migrate in these 
models of static disks. In fJSJ we discuss potential com- 
plication introduced by layer structure in the dead zones 
and feedback. In we discuss how this work connects 
to the ongoing effort of population synthesis modeling. 
In fj7]we summarize our results and discuss their impli- 
cations. 

2. MIGRATION RATES 

The type I migration rate is sensitive to the local E 
and T as well as to the radial gradients of these quan- 
tities. Following PBK10, we adopt the following general 
expression for the rate of type I migration, 



Reference 



Criteria 



dr . M p Er 2 



(2) 



where M p and M* are the planet and stellar mass, re- 
spectively, E, T, c s , and J7k all refer to the those quan- 
tities at the location of the planet. The magnitude 
and sign of the dimensionless efficiency factor / is a 
function of the dimensionless logarithmic surface den- 
sity and temperature gradients, p(r) = d log E/d log r 
and q(r) = d\ogT/d\ogr, respectively. Note that the 
PBK10 a and j3 parameters correspond to —p and — q, 
respectively, in this paper. 

In weakly perturbed disks, p and q are intrinsic prop- 
erties of the disk alone. The efficiency factor / is also a 
function of the viscous and thermal saturation parame- 
ters 



and 



p x = (n K r 2 x 3 s /27r X ) 1/2 = 3p,/2P r 1 / 2 



(3) 



(4) 



where Re = ^Ikt 2 '/^ttv is the Reynolds number, x s ~ 
(M p r /M^k) 1 / 2 is the dimensionless width of the horse- 
shoe orbit, h = Csfl^ 1 is the density scale height in the 
direction normal to the plane (i.e. the thickness) of the 
disk, and P r = v\ x is the Prandtl number. The magni- 
tude of both p v and p x depend on the mass ratio M p /M t 
as well as disk parameters, such that 



oc a-V 2 (AVM*) 3 / 4 (Vr)~ 7/4 



(5) 



in a disk in which the viscosity is parameterized accord- 
ing to eq. Q] 

Based on a 3D linear analysis of wave exc itation in 
a locally isothermal disk, iTanaka et al.l (j2002L hereafter 
TTW02) determined the migration rate to be 



fTTw(p, q) = 1.08(p + 0.8% - 2.52). 



(6) 



In typical disks, differential Lindblad resonances gen- 
erally lead to inward migration. Under some con- 
ditions, contribution from co-rotation resonances and 
horseshoe torques can le ad to outward migration. 
iPaardekooper et al.l (|2010al hereafer PBCK10) investi- 
gated disks in which there are fully- unsaturated (i.e. with 
a maximum efficiency) horseshoe torques. From these 



Radial Drift 
TTW02 

PCBK10 Locally Isothermal 
PCBK Adiabatic 



p > 1.5 - 0.5q 
p > 2.52 - O.Sq 
p > 0.85 - 0.9g 
p > 0.26 + 1.21q 



TABLE 1 

Conditions for Outward Migration 

upper torque limits, we can identify disk regions where 
planets have a potential to undergo outward migra- 
tion. In a similar (to TTW02), locally isothermal limit, 
PBCK10 found that the combined contribution from lin- 
ear Lindblad torques, linear-entropy related co-rotation 
torque, and non-linear vortensity related horseshoe drag 
(with a softening factor b/h = 0.4) leads to 



f iso (p,q) = -1.7 -1.8q + 2p. 



(7) 



PBCK10 also considered the effects of disk thermal 
structure. They found that in an inefficiently radiating 
disk with an adiabatic equation of state, a planet migrat- 
ing under the fully unsaturated (i.e. with full strength) 
non-linear horseshoe torques (with the adiabatic index 
7 = 1.4) yields 



fad(p,q) 



-1.22- 5.64g + 4.66p. 



(8) 



With these migration laws, it becomes simple to iden- 
tify the necessary conditions for stalling migration. The 
necessary criteria for outwards migration as a function 
of p and q are listed in table [T] For comparison and il- 
lustration purposes, we also included in this table the 
criteria relevant to the radial drift of pa rticles due to 
aerodynamic drag (jWeidenschillinglll977al ) . 

The results in table [1] indicate that with an adiabatic 
equation of state, outward migration in disks with neg- 
ative surface density gradients (p < 0) is only possi- 
ble if there is a strong negative temperature gradient 
(q < —0.22). Furthermore, in order for the torques to 
reach these fully unsaturated values, angular momentum 
must be extracted from the horseshoe region at the opti- 
mal rates, i.e. fast enough to avoid torque saturation, but 
slow enough to maintain a gradient across the horseshoe 
region. These requirements determine a set of sufficient 
conditions for stalling type I migration. 

3. DISK MODELS 

Here, we highlight migration rate's sensitive depen- 
dence on the disk structure. In the introduction, we have 
already indicated that the E and T distributions are de- 
termined by the not-fully-understood physics associated 
with the MRI turbulence and surface irradiation. In or- 
der to identify some robust features and assess range of 
uncertainties, we adopt three different disk models in- 
cluding 1) a relatively simp l e self - consistent constant a 
model (jShakura fc Sunvaevl [19731 : llir] fl98l (see 33~2j) 
and 2) a model for the irradiated outer regions of the 
disk (see M3.3p . We also consider 3) a disk model which 
takes into account the influence of a dead zone (see ^3.4|) . 

Before describing these three models, we orient the 
reader by describing some standard disk models and fea- 
tures common to most disk models. In figure Q] we show 
where various disk models lie on the p-q plane. The 
E profile of the minimum mass nebula corresponds to 
p = —1.5, and on this plot we place a star to indicate a 
geometrically flat (q = —1) MMN disk. It is interesting 
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Fig. 1. — The surface density (p) and temperature (q) gradi- 
ents for various disk models described in the text (Star-MMSN, 
Squares-CG97 , Circles- Viscous Disk, Triangles-GL07 . Diamonds- 
see i|3,4|l compared to the criteria for zero-migration under the 
different prescriptions presented in table [T] (Solid - TTW02, dotted 
- Aerodynamic Drag, dashed - Isothermal PBCK10, dash-dotted - 
Adiabatic PBCK10). The shaded rectangle indicates the location 
of constant a models. 



to note that in a MMN profile, the vortensity gradient 
is zero and the horseshoe toques do not contribute (see 
further discussion in Sj5j. 

The observed spectral energy distributions (SED) of 
young stars indicate that at least the outer regions of 
most protostellar disks have a flaring structure. This fea- 
ture means that h/r increases with r, an indication that 
these passive disks ha ve slowly declining tempe rature 
gradients (q > -1). (|Chiang fc Goldreich] 119971 here- 
after referred to as CG97) construct a model to match 
the observed features. In their model, there are three re- 
gions. Far from the star the disk is optically thin to both 
its own radiation and the stellar radiation, while in the 
inner region is optically think to both. In between these 
two boundaries, there is a marginally opaque region in 
which the disk is optically thick to the stellar photons, 
while optically thin to the longer wavelength disk pho- 
tons (see §3.3[) . The square symbols in Figure Q] indicate 
the temperature gradients of these three regions with an 
assumed MMN E profile. In this figure, the opaque and 
optically thin limits are virtually indistinguishable. 

3.1. Steady state disk models with constant a 

Instead of assuming a fixed E profile (such as that 
inferred empirically with the MMN scenario), it is also 
possible to deduce it from accretion disk models. In an 
viscous axisymmetric Keplerian disk, the mass accretion 
rate through each annulus is given by 



We assume short-period super-Earths were formed 
within a few AUs prior to their type I migration. In 
these inner disk regions, gas attains a quasi steady state, 
i.e. its mass accretion rate (M) is independent of disk 
radius (jRuden fc Linl |1986[ ) albeit M does evolve over 
time. In such a state, E is inversely proportional to the 
effective viscosity v. 

It is generally assume that turbulence provides an 
effective angular momentum transfer mechanism and 
the magnitude of viscosity can be approximated with 
an ad hoc a prescription in which v = ac s H where 
c s = (i? gas T//i) 1 / 2 (i? gas and fi are the gas con- 
stant and the mean molecular weight), H = c s /Qk 
and f^K = (GAfn./r 3 ) 1 / 2 are the sound speed, den- 
sity scale height and angular frequency of the disk 
gas (jShakura fc Sunvaevl Il973f ). If the efficiency fac- 
tor a is constant with radius, this steady state would 
lead to a simple relationship between the surface den- 
sity gradient and temperature gradient, such that 
p = —q — 3/2. The shaded rectangle in figure Q] indicates 
the parameter space occupied by steady-state, constant 
a models. 

In order to deduce a E distribution, it is necessary to 
independently or concurrently determine a T distribu- 
tion. One possible approach is to assume an equilibrium 
T oc r -1 / 2 distribution as i n the MMN model an d de- 
rive a E oc 7*" 1 distribution (Hartma nn et al.lll998| ). Al- 
though this equilibrium T distribution is adequate for the 
optically thin outer regions of the disk, it is not appro- 
priate for the inner disk regions which is mostly heated 
by viscous dissipation. In standard accretion disk mod- 
els, both T and E distribution can be obtained self con- 
sistently under the assumption of thermal equilibrium, 
i. e. viscous dissipation is balanced by local radiative flux 
(|Shakura fc Sunvaevl H973I : ILip|[l98lh . 

It is interesting to note that, although these standard 
disk models all lie in the lower region of this parameter 
space, the entire region indicated in this diagram satis- 
fies the stability criteria: d/dr[(r 2 £l) 2 ] > (Raylcigh's 
criterion) and dp/dr > dp/dz. Thus, the replacement 
of a constant a by a more general prescription does not 
generally introduce dynamically unrealizable models. 

3.2. Description of viscously-heated disk model 

In order to determine temperature of a viscously 
heated disk, we balance the radiative losses at the surface 
(of temperature T e g) with viscous heating and heating 
due to the background radiation, meaning 
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(10) 



where ctsb is the Stefan-Boltzman constant and Xj, is the 
background temperature. The effective temperature can 
be related to the midplane temperature via the effective 
optical depth T c 4 ff r G ff- = T 4 , where 



T e ff 



(11) 



M = 3twE 



37rac 2 f2 K 1 E. 



(9) 



3 V3 1 

-T+ 1 

8 4 4r 

and t = kE/2 is the true optical depth ()Hubenvlll990D . 
For all of the disk models we parameterize the opacity as 

k(T) - k T?. (12) 
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Fig. 2. — Evolution of transition radii as a function of mass- 
accretion rate. The transition between [i] viscously heated, [ii] 
strongly opaque, [iii] weakly opaque, [iv] marginally opaque, and 
[v] optically thin are marked with the dashed-dot, solid, dashed, 
and dotted lines. 

3.3. Description of GL07 model of irradiated disks 

In the determination of the E and T distribution for the 
intermediate and outer disk regions, contribution from 
stellar irradiation needs to be taken into account. In an 
attempt to introduce a more comprehensive prescription, 
iGaraud fe Linl (|2007| ) constructed an analytic disk model 
(hereafter referred to as GL07) based on the passive disk 
model of CG97. However, instead of assuming a MMN 
surface density profile, GL07 calculated the structure of 
a disk evolving with a constant a-viscosity. Addition- 
ally, this model includes both stellar irradiation and vis- 
cous heating. This prescription results in five distinct 
regions. The outer three regions, the optically thin [v], 
marginally opaque [iv] , and weakly optically thick [iii] re- 
gions are similar to those in CG97. In these regions the 
T distribution, is independent of the E profile, although 
a steady state solution can be preserved with a variable 
a prescription. 

In the GL07 model, there are two regions close to the 
star , that do not occur in the CG97 passive-disk model, 
a strongly-opaque region [ii] and viscously-heated region 
[i]. In region [ii], the optical depth to the disks own ra- 
diation is high, causing the disk to be non-isothermal 
both in the disk plane and in the direction normal to it. 
In region [i], viscous dissipation dominates the thermal 
profile. In these regions q is a function of p. In the in- 
ner disk regions where accretion flow is in a quasi-steady 
state, both quantities (p and q) are well defined. 

These five regions ( [v] and [iii] are nearly indistinguish- 
able) are indicated by triangles in figure [T] The radial 
transitions between these regions evolve as a function of 
the mass accretion rate (see figure [5]) . It is interesting 
to note that the optically thick region of GL07 lies re- 
markably close to the separatrix between outwards and 
inwards migration in figure [TJ 

3.4. Description of an MRI-active Disk Model 

As we have indicated above, the T distribution in the 
inner disk regions is determined by viscous dissipation 
and that in the outer disk regions is determined by stel- 
lar irradiation. The evolution of the E distribution ev- 



erywhere in the disk is determined by the efficiency of 
angular momentum transfer, i.e. the magnitude of ef- 
fective viscosity v, which we assume is mainly due to 
MRI driven MHD turbulence. In the innermost regions 
of protostellar disks where hot Jupitcrs and super Earths 
eventually reside, the ionization fraction of the disk gas is 
modest and MHD turbulence is effective throughout the 
disk. But, beyond a fraction of an AU where most planets 
form and start their migration, disks' midplane is gener- 
ally too cool and neutral to be affected by magnetic fields. 
However, their surface layers are exposed to the stel- 
lar UV photons and sufficiently ionized for MRI-driven 
MHD turbulence to prevail. Wave propagation from the 
surface layer to the midplane can also le ad to angular mo- 
mentum, albeit at much reduced rates ((Fleming fc Stond 
l2(ffl!lTurner et al.ll2n07Tl . 

Dynamical evolution of MHD turbulence in protostel- 
lar disks is beyond the scope of this paper. In lieu of a 
detailed turbulent-disk model, we simply adopt a vari- 
able a in disk regions extend over the migration paths 
of hot Jupiters and super Earth. This ad hoc proce- 
dure is const ructed in the same spirit as the original a 
prescription (jShakura fc Sunvaevlll973T ). i.e. for compu- 
tational simplicity. To do so, we introduce a parame- 
ter £(r) = dlog a/d log r which varies with r depending 
on the ionization of the disk gas. In the appendix, we 
present a simple modification to the GL07 model to allow 
for a variable a. 

In locations where thermal ionization becomes suffi- 
cient to sustain MRI (~ 10 3 K) and at the sublimation 
front where dust destruction decreases the recombina- 
tion rate, the magnitude of a increases rapidly with the 
disk temperature. If viscous dissipation is the dominant 
contributer to the temperature profile in these locations 
(e.g. in the regions [i] and [ii]), the disk would tend 
to evolve and attain a relatively shallow radial T gra- 
dient. In the limit that q is negligible, (y = —9/2 and 
Cso = ~ 3, leading to a py = 3 and pso =3/2 for the vis- 
cous and strongly opaque regions, respectively. In both 
cases, fully unsaturated horse shoe torque would induce 
a rapid outward migration. 

3.5. Inner Disk Truncation 

A generic feature of protostellar disks is the presence 
of a cavity at their center which is cleared by the mag- 
netosphcre of their central stars. The presence of this 
magnctosphcric cavity quenches both Lindblad and horse 
shoe torques and therefore provides a natural post fo r 
hot Jupiters' type II migration to stall (|Lin et al.lll996| ). 
Here, we consider the possibility that the termination 
of super Earths' type I migration may have also been 
strongly influenced by this structure. 

As an additional feature in our disk models, we con- 
sider the possibility of an inner boundary due to the 
disks' interaction with the stellar magnetosphere of their 
central stars. The location of this magnetospheric trun- 
cation radius is assumed to occur at 

r mag ~ 0.5^ /7 (2GM st ) 1 / 7 Af- 2 / 7 , (13) 

wher e /x* is the stellar magnetic dipole moment (|Konigll 
fl99l . 

Many of the approximations integral to the disk model 
break down in this location. The actual E and T profiles 
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in this disk region need a detailed analysis of the process 
of mass loading onto the stellar magnetic field as well as 
magnetic diffusivity and dissipation near the disk inner 
boundary. These studies are beyond the scope of this 
paper. For the purpose of studying the destiny of super 
Earths, we impose a sharp positive surface density gra- 
dient near r mag , where the surface density is reduced by 
the factor / mag - In this paper we use 

Hte) 2 )' (i4 » 

though we n ote that this is an arbitrary choice. And 
as noted by jTsang 20Tll) the interaction at the inner 
edge of the disk are likely more complicated than the 
simple picture presented here. As a planet approaches 
the inner edge of the disk, it is possible that reflections at 
the inner disk edge may cause the planets to have their 
migration halted at a few times the inner disk radius. 
But while the details arc different, the unsaturated horse 
shoe torque near r mag generally induce a rapid outward 
type I migration, creating a planet trap. 

4. TYPE I MIGRATION RATES DEDUCED FOR VARIOUS 
DISK MODELS 

4.1. The Migration Rate of a Planet in a Steady-State 

Disk 

In order to discuss the dominant and relevant physi- 
cal effects, we illustrate the migration of a planet in a 
static (i.e. a steady-state) disk. In figure [3j we show 
the relative migration rates (as compared to the TTW02 
model) for a disk which is heated entirely by viscous dis- 
sipation (see ^3.2|) . Such a structure is possible if the 
stellar irradiation is bl ocked by self shadowing effects 
(jDullemond et al.ll2002f h We adopt the following model 
parameters: M = 10- 8 M Q yr- 1 , a = 10" 3 , T b = 20K, 
and r mag = 0.15 AU. 

From this standard viscous disk model, we find dis- 
tinct regions where the disk structure can allow outwards 
migration. Other than the inner boundary of the disk, 
these locations include the regions far from the subli- 
mation temperatures, i.e. zone 1 where the opacity is 
dominated by silicates (f3 = 1/2 in eq IT2"j) and zone 2 
where the opacity is dominated by ice grains (f3 = 2). 
In the regions interior to zone 1 (where silicate grains 
partially sublimate), between zones 1 and 2 (where ice 
grains partially sublimate), and outside zone 2 (where 
the disk become optically thin in the direction normal 
to the disk) , fully unsaturated horse shoe torques cannot 
offset the differential Lindblad torque to induce planets 
to migrate inward. . 

4.2. Potential location of planet traps 

Fully unsaturated horse shoe torque has a tendency to 
induce planets formed in zones 1 and 2 to converge with 
those formed beyond these zones. However, due to satu- 
ration effects, only planets with a limited range of masses 
can experience the full horse shoe torques. For example, 
the dynamics of very low-mass planets is not affected by 
the horse shoe torques and they migrate inwards every- 
where in the disk, albeit slowly, under the influence of 
differential Lindblad torque. Horse shoe torques are also 
saturated for relatively massive planets (with masses at 



least an order of magnitude larger than that of Earth) 
and they undergo rapid inward migration everywhere ex- 
cept near the disks' inner boundary. Their inward mi- 
gration would be stalled at r mag if the surface density fall 
off rapidly interior to the edge of the disks' inner cavity. 

4.3. Migration during planetary growth 

We now consider the implication of the results in fig- 
ure [3] on the growth of planetesimals and protoplanetary 
embryos. For this model (i.e. this choice of M and a), 
all planetesimals with mass larger than that of Mars mi- 
grate inward at rates comparable to that obtained by 
TTW02. However, with adequate metallicity, the charac- 
teristic growt : htime_scale for the modest-mass planetes- 
imals (|Ida fc Linl I2004T ) is shorter than their migration 
time scale such that they can grow in situ into dynami- 
cally isolated protoplanetary embryos. 

Embryos' isolation mass increases rapidly with r. Out- 
side an AU or so, the it exceeds Mars' mass. Based on 
the results in figure [H we find that the partially satu- 
rated horse shoe torque induced on these modest-mass 
embryos is adequate to cause an outward migration to- 
wards the outer boundaries of zones 1 and 2. 

The migration of these embryos would be stalled at 
these planet-trapping boundaries. In this set of mod- 
els, there is an arbitrary growth barrier. Once a core is 
stalled at a barrier, it would accrete all of the material 
in its vicinity until it has attained a new isolation mass 
(Mj so is modified by the enhanced £^ due to the accu- 
mulation of planetesimals and embryos) . In a disk with a 
single core there is no physical mechanism for it to grow 
beyond this size. 

In this model where the disk is entirely heated by vis- 
cous dissipation, it is possible to trap super Earth em- 
bryos beyond a few AU. Such embr yos may have masses 
comparable to that (M c ~ 10M®, [Pollack et all [19960 
required to initiate efficient accretion of gas onto the pro- 
genitor cores of proto gas giant planets. 

The actual mass range and radial extent of zones 1 and 
2 depend on the magnitude of M. During the depletion 
of the disk, zones 1 and 2 cover smaller planetary masses 
and radial distance from their host stars. 

4.4. Sensitive dependence on the disk models 

In figure [31 we show the landscape for a GL07 disk 
with M = 10- 8 M o yr- 1 , a = 1CT 3 , = 1. In the in- 
ner most region of the disk viscous heating dominates, 
and it is analogous to zone 1 in the fully viscous disk 
model (see discussions in the previous section). Conse- 
quently, planets with masses in the appropriate range 
undergo outward migration. Outside this region, stellar 
irradiation provides most of the heating for the disk and 
modifies the temperature distribution. Consequently, the 
horseshoe torque associated with the co-rotation reso- 
nance is weakened. There is a distinct strongly-opaque 
region where the horse shoe torques reduce the impact 
of differential Lindblad resonance and significantly slow 
down but does not reverse planets' inward migration, 
even for those with optimum masses. The main reason 
for this reduction in the inward migration rate is due to 
the relatively steep decline of the mid plane temperature 
with the radius. With the inner region being heated by 
viscous dissipation, the temperature gradient i.e. q(r) in 
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Fig. 3. — The temperature (top), surface density (middle), and 
migration rate relative to the magnitude of TTW02 as a function 
of radius and planetary mass in a viscously heated disk (bottom). 
We hatch the region in which migration is outwards for emphasis. 
The magnetospheric truncation radius is indicated by the dotted 
vertical line. The black dashed-li nc i ndicates when the disk in- 
teraction becomes non-linear (eq. 1161) and the white dashed-linc 
indic ate s the planet mass which satisfies the gap opening criteria 
(cq. I17p . In this example the two criteria are virtually identical, 
but this is not the case in general disks. 

GL07 is more negative than that for a totally passive 
passive, flaring disk in CG97. 

For these choices of parameters, the most massive em- 
bryos that can be trapped is ~ 3 — 4Mq in an unper- 
turbed disk at a fraction of an AU. This mass range is 
well below that (M c ) required for the acquisition of mas- 
sive gaseous envelope for gas giants through the onset 
of efficient gas accretion. The location is also too close 
host stars to be compatible with that of most known gas 
giants. The upper mass limit of retainable embryos in- 
creases with the poorly determined and constrained mag- 
nitude of a. But, for a given M, the surface density of the 
gas and presumable the heavy elements decreases with a. 
A larger value of M or a smaller value of a may also lead 
to the growth and retention of more massive cores. The 
torque applied by the cores can also lead to modifications 
in the £ distribution and this feed-back effect can also 
reduce the type I migration rate (see $5]). 

4.5. Disks with dead zones and active layers 

In figure [5] we show the relative migration rates for a 
planet in a disk with a MRI active surface layer. For 
simplicity, we first consider a model in which the surface 
density of the MRI active layer is a constant (XU = 
50gcm -2 ), similar to the seminal [Gammiel (|1996l) paper, 
and potentially there is a dead zone below it with surface 
density Ed . We take the viscosity in the active layers to 
be large with an a a = 10 ~ 2 , while the viscosity in the 
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F ig. 4 . — The same as figurc|3]for the GL07 model, as described 
in jOl 

dead zone is small old = 10~ 4 . With this small but non- 
vanishing value for old, we approximate the propagation 
of turbulent disturbances near the disk surface into the 
midplane regions of the disk. The ac tual values of 
and an are still poorly constrained (jFleming k, Stonel 
l200llTurner et al J 120071) . 

In regions where £a is smaller than the total surface 
density, accretion flow in the active surface layer of the 
disk leads to a smaller effective magnitude for the verti- 
cally averaged a = {ola^a + £*£)££))/£. The main con- 
sequence of this modification is the expansion in both 
the lower and upper mass range for embryos that can 
be retained. The upper mass range exceeds the magni- 
tude of M c . The lower mass range extends well below 
the isolation mass for embryos even at relatively small 
radius. However, if we take into account the effect of 
vertical structure on the saturation parameter, the lower 
mass range would be limited to a fraction of the Earth 
mass. Although these tendencies maybe favorable to the 
formation of gas giants, the location of the planet trap- 
ping radius is only a fraction of an AU. Larger values 
M are needed for the planet trapping radius to become 
comparable with the typical orbital semi major axis of 
known gas giants. 

5. SOME POTENTIAL COMPLICATIONS 

A robust feature in the above results is a tendency for 
super-Earths to migrate outwards in the inner regions 
and inward in the outer regions of the disk. This property 
implies the presence of a planet-trapping radius where 
these planets may converge and accumulate. We also 
showed that the mass range of the potentially trapped 
planets and the location where their type I migration is 
stalled are model-dependent and subject to uncertainties. 
In this section, we list several additional subtle effects 
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dr/dt relative to |TTW| dr/dt relative to |TTW| 



Fig. 5. — The same as figure Q5] for a disk with a dead zone and 
an active layer with a constant column depth (see t|4.5l) . 

which may also affect these quantities. 

5.1. Disks with variable active layers 

In figure [6l we show the relative migration rates for a 
planet in such a disk with a T,a which varies with the 
local dust -to-gas ratio. We use a n approximate fit to the 
results of lllgner fe Nelsonl ([2006D [model 4] in which 

where fd is the dust-to-gas ratio for micron-sized par- 
ticles. This prescription introduces a feature at the 
snow line where the dust to gas ratio decreases as a 
planet moves inward s due t o the sublimation of ice grains 
(jKretke fe Lml 120071 120 lOf ). In this confined transition 
region, embryos can be trapped at several AUs from 
their host stars. However, we note that for these pa- 
rameters the mass of cores that can be trapped is rela- 
tively small, below M c ~ 10M®. Whether or not these 
cores are sufficiently massive to accrete atmospheres and 
become gas giants depends upon properties such as the 
pla netesimal mass accret i on rate and a tmospheric opac- 
ity (IHubickvi et al.ll200l IRafikovi[2006l) . 

5.2. Revisiting the Horseshoe torques in the Dead Zone 

In addition to the reduction in a and the correspond- 
ing modification in the £ and T distribution, a layered 
disk structure may also affect the saturation of embed- 
ded planets' co-rotation torque. The horseshoe region 
where the co-rotation torque is applied has a dimcn- 
sionless width in the disk plane x s ~ (Mpr/M^h) 1 ^ 2 
(|Paardekooper fc Papaloizoul2009t) . In the direction nor- 
mal to the disk plane, the horseshoe region extend to a 
similar dimcnsionless height z s ~ x s . 



Fig. 6. — The same as figure [3]for a disk with a dead zo ne w ith 
variable thickness depending on the dust to gas ratio (sec ^5. It . 

For a 10 M e planet (M p /M Q ~ 3 x 10~ 5 ) embedded 
in the midplane of a disk with a typical h/r ~ 0.05, z s ~ 
0.02. For these parameters, there is a substantial dead 
zone that extends over one disk scale height (so that the 
dimcnsionless height of the dead zone Zd > h/r) and the 
entire horseshoe region is engulfed within the dead zone. 
In this situation, the horseshoe region is entirely isolated 
from the active layer. Ideally, the continuous viscous 
diffusion of disk gas into the horseshoe region of relatively 
low mass planets may be quenched if their z s < Zd- It is 
unclear if the Reynold's stresses which provide a source 
of angular momentum transport in the dead zone (and 
therefore generate an effective an) may also redistribute 
vortensity and entropy between the horseshoe region and 
the surrounding disk with a similar efficiency. Without 
such redistribution, the co-rotation torque of low-mass 
planets with coplanar orbits would be saturated after a 
few libration periods. 

In principle, the limited saturation for planets with 
h/r > z s > Zd can be estimated from 3D hydrodynamic 
simulations. But, such simulations are time consuming, 
especially for such low masses. Therefore we utilize our 
approximation to highlight how including this effect may 
impact the migration rates. In figure [7] we demonstrate 
the results if we assume that the effective viscosity is 
used in determining the Reynolds number (Re) in eq.[3]is 
negligible if z s < Zd, meaning that, for the sole purpose of 
calculating Re, we set a = 10~ 8 . Under this assumption, 
the horseshoe torque saturate for nearly the entire dead 
zone, resulting in inwards migration for almost all masses 
and positions. 

As a comparison, in figure [5] we assume that if z s < 
Zd then the relevant effective viscosity is the dead zone 
viscosity ar> ■ In this case the horseshoe torques saturate 
at a much smaller mass so we still see stalling locations 
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Fig. 7. — The same disk model as shown in figure [5] but with the 
assumption that the co-r otation torques fully saturate in the dead 
zone if z a < z^ (see SI5.2I1 . 
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Fig. 8. — The same disk model as shown in figure \E\ for the 
assumption that the effective viscosi ty to desaturate the corotation 
torques is equal old if z a < z^ (sec ij5.2t . 

for small planctcsimals, but we lose the ability to stall 
cores greater than around 3 M®. 

In these estimates, we assume the cores have co-planar 
orbits. In N-body simulation of planetcsimal dynamics 



and embryos formation, their mutual gravitational per- 
turbation and collisions can lead to dispersions in the 
velocity (er z ) normal to the disk plane. During the initial 
growth, <T Z of the most massive embryos is damped by 
dynamical friction. But during the advance Oligarchic 
growth stage, perturbation between comparable mass 
embryos may induce excitation o z to be a fraction of their 
surfa ce escape speed, despite gas drag (|Kokubo fc Idal 
1998). If o z > f^K^d, the embryos would encounter MRI 
active regions during at least a fraction of each Keplerian 
orbital period and the reduction of the diffusion through 
the horseshoe regions would no longer be an issue. 

5.3. Feedback and transition to type II migration 

In the above sections, we assume that planets with 
masses up to a few M ffi do not significantly alter the 
average disk E and T distribution. In our evaluation on 
the direction and rate of such planets' type I migration 
and the extent of saturation (i.e. the magnitude and sign 
of /), we calculate the E and T distribution neglecting 
any feedback effects. However, we note that previous 
research has demonstrated that, even with a modest mass 
(such as ~ 5M ffl ), planets may alter the disk structure, 
especially in the horseshoe region, an d slow down the 
rate of migration dMasset et al.|[2006al) . 

Additionally, the type I migration rate is modified once 
the interaction between the disk and a planet become 
non-linear. This limit is reached when the planet's hill 
sphere becomes comparable to the disk scale height, i.e., 
once 



Mp 



> 1. 



(16) 



Additionally the plan et begins to open a gap 
(|Lin fc Papaloi zou 1993]), when 



e 2 



40a 




< 1. 



(17) 



Gap formation interrupts accretion flow through the disk 
and conseque ntly triggers a transition from type I to type 
II migration (|Wardlll997D . This limit is more relevant to 
gas giants rather than super-Earths. 

Though super-Earths do not generally have sufficient 
mass to open gaps and, their feedback can significantly 
affect type I migration process. In the low viscosity (with 
a << 1) limit, gas in the horseshoe region encounters 
weak shock when it flows close to a-few- Earth-mass plan- 
ets' Roche lobe. Wave dissipation across the shock leads 
to angular momentum deposition an d evolution in t he E 
distribution in the horseshoe region (iLi et al.H2009h . On 
the horseshoe libration time scale, this effect can signif- 
icantly reduce t he torque magn itude and planets' type I 
migration rate (|Yu et al.ll20l"0ri . 

Shock induced dissipation along the stream lines also 
modifies the vortensity (the vorticity divided by the sur- 
face density) along each streamline. The development of 
extrcma in the vortensity distribution across the stream 
lines promotes the growth of secondary or Rossby Vortex 
instabilities (RVI) which can be suppressed by adequate 
viscosity. In the limit of negligible viscosity, the growth of 
RVI c an lead to the emergence of large vortices (|Li et al.l 
I2005T ) . The torque induced by the interaction between 
the planets and the vortices is chaotic. It can affect both 
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the direction and rate of type I migration. Nevertheless, 
the onset of RVI significantly re duces the effect ive mag- 
nitude of type I migration rate (lYu et alJl20iH . These 
feedback effects can extend the upper mass limit for the 
trapped planets. 

6. RELATIONSHIP TO POPULATION SYNTHESIS MODELS 

In order to assess how these various disk models may 
impact the emergence and appearance of planetary sys- 
tem one may wish to include these prescriptions in pop- 
ulation synthesis models. We note that without actually 
preforming a detailed population synthesis model we can- 
not compare how precisely the models will be affected, 
but we c an compare these m odel s to existing form ula- 
tions. In llda fe Linl (|2008aD and llda fe Linl pOloh the 
authors assume that the gas follows the radial surface 
density gradient of the minimum mass solar nebula and 
the temperature follows a power- law of T cx r -1 / 2 . As 
indicated in figure [TJ this disk structure would lead to 
inwards migration under both the isothermal and adia- 
batic regimes. Therefore if one assumes that disk model, 
including the additional thermal effect and saturation ef- 
fects discussed in this paper will not significantly change 
the result, that type I migration at its predicted values 
is too fast to form the observed planet population. In 
llda fc Linl (|2008bl ) where the authors adjust the surface 
density to include a simple prescription for the dead zone. 
However the same thermal profile is still used. In this 
case fully considering the possibly of torque saturation 
causes the snow-line "planet trap" to be only effective 
for a limited range of masses. The largest change for 
these models is that if they used a self-consistent tem- 
perature model they could create inner, viscously heated 
region. In the papers as given the authors recognize that 
in this region they must artificially stop the planets mi- 
gration, and the viscously heated inner region may be a 
natu ral way to achiev e this effect . 

In iMordasini et al.l (|2011[ ) th e authors u s ed th e vis- 
cous heating presc ription of iLvra et all (|2010l ) and 
ISandor et all (|2011h which is very similar to the model 
described in ^3.21 and figure [3] This leads to various 
"concentration points" for planet migration. However, 
in this paper we argue that a viscously heated model 
may not be valid for large regions of the disk. If indeed 
the disk is passively heated from the star then the disks 
will not have planet traps far from the parent star un- 
less there is dead zone or other change in viscosity. If 
this is the case than the results may b e more similar to 
those shown in IMordasini et ail ([2012D in which the au- 
thors use a viscously heated model but do not account 
for the thermal effects. This means that they see inwards 
migration throughout the disk without stalling zones. 

7. SUMMARY AND DISCUSSION 

In this paper we discuss the migration of low-mass 
planets (super-Earths) in gaseous disks via type I migra- 
tion inclu ding the impact of h o rsesho e torques as pre- 
sented by iPaardekooper et all ()2010bl ). We find that 
planets with a few Af^'s have a tendency to undergo 
outward migration through regions which are predomi- 
nantly heated by viscous dissipation and inward through 
regions which are mainly heated by stellar irradiation. 
Since typical protostcllar disks are mostly heated by vis- 
cous dissipation at small radii and by stellar irradiation 



at large radii, these planets tend to converge and stall 
at a fraction of 1 AU from their host stars. This result 
is robust and general for self-consistent protostellar disk 
structures in accordance with the conventional a mod- 
els. Thus, we anticipate a high retention efficiency for 
super-Earths. This extrapolation is consistent with the 
observed ubiquity of super-Earths. 

However, the location of the stalling radius is some- 
what model-dependent. Our numerical results indicate 
that its value is smaller in models that take into consider- 
ation stellar irradiation than in model s that are primar - 
ily heated by viscous dissipation (e.g. lLvra et al.|[2010t) . 
In the composite (viscous dissipation plus surface irra- 
diation) models, it is difficult to sustain outward mi- 
gration to large radii without any self-shadowed effects. 
Whi le self-shadowin g by a puffed up inner rim is possi- 
ble (|Ke et al.ll20ilh . it is unclear if this shadowing can 
extend beyond a few AUs. 

In all of the models the planets migrate inwards or out- 
wards until they reach a barrier such as inner-edge of the 
dead zone ry, the magnetospheric truncation radius, or 
other stalling radii. However, even if a planet stalls at a 
given radius within the disk it may not determine its final 
location because, as a disk evolves, the location of these 
barriers tend to change. For example, r mag increases and 

rid decreases with diminishing M. This sensitive depen- 
dence on the disk structure may lead to large dispersion 
in the asymptotic period of super-Earths. In contrast 
to a 3-day pile up in hot Jupiters' period distribution, 
super-Earth candidates, which are identified by Kepler's 
mission, show a nearly uniform logarithmic period distri- 
bution without any preferred orbital configuration. The 
observed spread-out period distribution of super-Earths 
does not contradict the expected large dispersion in both 
the disk structure and planets' stalling condition. 

For planets that do reach the inner-edge of the disk, 
the stalling locations for low-mass planets differs from 
those of more massive (Jovian mass) gap-opening plan- 
ets. These high-mass planets terminate their type II 
migrating once they have reached a 2:1 resonance in- 
terior to the inner edge of the truncated disk or in the 
limit t hat the disk ma ss becomes smaller than the planet 
mass (|Lin et al.l 11 9961 ). Therefore we expect that short 
period high mass planets (hot Jupiters), the product of 
type II migration, should be found on average interior 
lower mass planets, whose migration is stalled at the 
magnetospheric truncation radius. As we continue to 
detect low-mass planets on short-period orbits they will 
become a important tool to constrain models of planet 
formation and migration. The final stopping location 
of short-period planets will reveal information about the 
structure of the inner disk region. 

Additionally, as type I migration is proportional to the 
planetary mass, planets on the upper range of the type 
I migration regime are the most vulnerable. For these 
higher-mass embryos the horseshoe torque begin to sat- 
urate, meaning that they migrate inwards at the full rate. 
However, these embryos are also the planets which are 
able to modify the £ profiles of the gas disk in their 
vicinities. These feedback effect could quench type I mi- 
gration, induce vortex formation through RVI, or lead 
to type II migration. These effect may lead to the fall 
off in the observed size distribution (at around 3-4 Earth 
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radii). 

High mass-accretion rate promotes core retention at 
large radii. This correlation is primarily due to the 
fact that the aspect ratio (h/r) generally increases with 
M. In this limit, the horseshoe torque does not satu- 
rate until the embryos have acquired a relatively large 
planet mass. For sufficiently large M(> 10 _7 il/ Q yr _1 ). 
stalled embryos can attained sufficiently large mass to 
initiate efficient gas accretion. A larger a can also pre- 
vent inwards migration of these particularly vulnerable 
cores. However, the location of the transition from the 
viscously heated region to the strongly opaque region 
(i.e. the transition from inwards to outwards migration) 

depends only on M not on a. In general the existence of 
a dead zone hinders the formation of cores at large radii 
because the low viscosity prevents the barrier for high- 
mass cores. However, the additional degree of freedom 
provided by the variable a parameter open the possibil- 
ity of rapid outward migration for embryos with various 
masses. This effect may create planet traps, but only if 
the horseshoe torques can remain unsaturated. 

A model in which planetary cores are retained at the 
"planet traps" begs one very important question. Once a 



planetary embryo acquires more than one Earth mass, it 
would migrate to one of these trapping locations on rel- 
atively short timcscales, unless there is another mecha- 
nism to slow type I migration. Therefore, if, as expected, 
a system of multiple embryos emerge in a disk around a 
common host star, they would have a tendency to mi- 
grate to locations which are already occupied by others. 
This convergent dynamical evolution raises the question 
of how two (or more) embryos interact as they approach 
one another. This interaction needs to be explored in 
future studies, especially as it is clear that super-Earths 
appear to primarily be found in multi-planet systems. 

The prescription generated in t his paper can be applie d 
to population synthesis models (jlda fc Linl 12004 12010ft . 
Through statistical comparisons between these models 
and the observed data, it is possible to highlight both 
robust physical processes which lead to common features 
and marginal effects which introduce diversity and dis- 
persion. 
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APPENDIX 
GENERALIZATION OF GL07 

While GL07 is presented for a disk model with a constant a parameter, it can be generalized to include radial 
dependence. The temperature models presented in GL07 are entirely local, and therefore independent of the radial 
variation of a. So for an arbitrary surface density profile one may simply use the prescriptions presented in GL07. 
However, in order to determine the surface density profile, assuming a steady-state disk, 

P = -q-C~l (Al) 

We define £ = din a/ In r. In the optically thin, marginally opaque, and weakly opaque regions the temperature 
profile is independent of the surface density profile so one can directly compute the surface density profile from this 
relationship. In the strongly opaque and viscously heated region the temperature gradient depends on the surface 
density 

15n + 15 - 9/3 + C(14n + 14 - 2/3) , 
PSO = 2(7n + 15-/3) (A2) 



where n is the polytropic index and the opacity varies according to k(T) = Ky(T/T*) 



9n - 3/3 + 15 - 4C 

iv = i^Ts 3 (A3) 



For the fiducial model in which n = 2 and ft = 1, 



and 



In the viscously heated region 



and 



PSO - -!±££, (A4) 
,so - -S±£. (A5) 
Pv = t, — (A6) 



(AT) 



We note that these are valid only under the same limit of GL07 that 



r amr n z ar \ r 
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